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The anomalous processes, 7 — > Stt and 771 tttt, are investigated within the Dyson-Schwinger 
framework using the rainbow-ladder approximation. Calculations reveal that a complete set of 
ladder diagrams beyond the impulse approximation are necessary to reproduce the fundamental 
low-energy theorem for the anomalous form factor. Higher momentum calculations also agree with 
the limited form factor data and exhibit the same resonance behavior as the phenomenological vector 
meson dominance model. 
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INTRODUCTION 



Anomalies, in particular the anomalous nonconserva- 
tion of the axial vector current, are a consequence of 
renormalization and are thus potentially present in all 
gauge theories. In fact a conserved axial current is in- 
compatible with electromagnetic gauge invariance since 
radiative corrections lead to a non-zero current diver- 
gence. The quintessential example is the 7r° 77 de- 
cay, which, in the chiral limit, would be forbidden if the 
axial current was conserved. However, the form fac- 
tor for this process, F'^'^iQ'^), is non-zero and in the 
combined chiral and zero-momentum limit it is given 
by ^ ^2^(0) = eV(47r2A) = a/{nU) 

where e is the proton charge, a is the fine structure con- 
stant and /tt is the pion decay constant in the chiral limit. 

This work addresses another anomalous process, 7 — > 
tt+tt'^tt^, which also should, but doesn't, vanish in the 
combined chiral and soft momentum limits. The form 
factor for this, and the crossing related process 777 —>■ tttt, 
is denoted by i^'^'^(s, u), where s = —{Ps + Pi)^, 
t = -{P2 + Pa?, and u = -{P2 + P3? . with P, 
the pion momenta. In the chiral limit, and for zero 
four-momenta, Refs. HSU proved #3^(0,0,0) = 
F^V{efl) = e/{A^^fl). 

Several issues, both experimental and theoretical, mo- 
tivate this investigation. Experimentally, only limited 
low momentum information about F'^'^ is available from 
Primakov production using the ttA — ^ tttt A reaction 
and, for timelike photons, from e^e~ — > Stt measure- 
ments IM 0- However, measurements to determine 
the form factor for the process 77r — s- tttt are now under- 
way at JLab . New results for the form factor F^"^ in 
the range 0.27 GeV^ < s < 0.72 GeV^ are expected to 
be released in the near future pd| . 

Theoretically, there are published analyses [T^ IT^ IT^ 
based on the set of Dyson-Schwinger equations [DSEs] in 
the rainbow-ladder truncation. These pioneering studies 
evaluated diagrams in the generalized impulse approxi- 



mation [GIA]. However, it is necessary to go beyond the 
impulse approximation in order to correctly describe tt-tt 
scattering : consistency requires inclusion of all 

planar diagrams ignoring gluon self-interactions (triple 
and quartic gluon vertices). In the current study we ap- 
ply the same scheme to 7-877 processes and elucidate that 
the rainbow-ladder truncation in combination with the 
GIA is insufficient for a realistic description. 

The structure of this paper spans six sections. In the 
next section the DSE formalism for mesons is applied to 
evaluate F^"^ . Our numerical results for the form factor 
for symmetric pion momenta are presented in Sec. IIIII 
and the low-energy theorem is reproduced. In Sec. IIVI 
we compare our results with meson exchange models for 
more general pion momenta. Finally, we confront the 
limited data in Sec. and then summarize our conclu- 
sions in last section. 



II. MESON INTERACTIONS IN THE 
DYSON-SCHWINGER APPROACH 

A. Dyson— Schwinger equations 

Here we briefly summarize the DSE formalism applied 
to mesons. For more complete details consult Refs. 0, 
0, 0, li^. Working in Euclidean space, with met- 
ric {7^, 7j.} = 25 fj,^, jI = 7p and a-b = a^bi = X^Li O'l^i^ 
the DSE for the renormalized quark propagator having 
four-momentum p is 

Z, / g^D^,ik)^j^S{q)T^{q,p). (1) 
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Here a — 1...8 is the color index, Dfj,^{k) the 
dressed-gluon propagator and F" {q, p) the dressed-quark- 

gluon vertex. The symbol represents ^.^^^4 
with k — p ~ q. Equation ^ has the general solu- 
tion S{p)^^ — i^A{p'^)+ B{p^) which is renormalized at 
spacelike so that A{^^) = 1 and B{fi^) = mq{ii) with 
niq^fi) being the current quark mass. 

Mesons are modeled as bound states governed by 
the Bethe-Salpeter amplitude [BSA], Th, which is a solu- 



tion of the homogeneous Bethe-Salpeter equation [BSE] . 
For quark flavors a, b and incoming, outgoing momenta 
p+=p + P/2, p- = p- P/2, the BSE is 

r'S{p+,P-) = rKip,q;P) 

Jq 

xS'^{q+)r%\q+,q_)S\q^). (2) 

The momenta q± are similarly defined. The renormal- 
ized, amputated qq scattering kernel, K, is irreducible 
with respect to a pair of qq lines. Equation ^ only ad- 
mits solutions for discrete values of = —m^, where 
■niH is the meson mass. Imposing the canonical normal- 
ization condition for qq bound states then uniquely de- 
termines Tff. For pseudoscalar bound states the most 
general decomposition for Tps = is plll2^ 

rps{q+,q-)^"f5[^Eiq^q-P) + rFiq^,q-P) 
+ iG{q\q-P)+a,,P^q'' H{q\q-P)\ , (3) 

where the amplitudes E, F, G and H are Lorentz scalar 
functions of and q ■ P. The odd C parity of the neutral 
pion requires that the amplitude G be odd in q- P, while 
the others are even. 



B. Rainbow- ladder truncation 

The central approximation in this work is the rainbow- 
ladder truncation for the set of DSEs. For the quark 
DSE, Eq. the rainbow truncation is 

Zig'D,Uk)K{q,P) g{k')D'^^f{k)j,^, (4) 

where D^^f,'^{k — p ^ q) is the free gluon propagator in 
the Landau gauge and G{k'^) is an effective qq interaction 
that reduces to the perturbative QCD running coupling 
in the ultraviolet region. For the BSE, Eq. the ladder 
truncation is 

K{p, q- P) ~g{e) D'p^k)^^,^j, . (5) 

These two truncations, in tandem, yield vector and 
axial vector vertices satisfying their respective Ward- 
Takahashi identities. This ensures chiral symmetry is 
respected with the pion being a massless Goldstone bo- 
son in the chiral limit ^2ii i22e 1^ 113 - Further, a con- 
served electromagnetic current is also obtained if the GIA 
is used to calculate electromagnetic form factors . 

C. Amplitude for "/-3n processes 

The above formalism is now applied to the 7 ^ Stt 
process to compute the anomalous amplitude, A^'^ , and 
the corresponding form factor, F^'". For simplicity, let 
us first consider one of the configurations that contribute 
to this process, depicted in the GIA in Fig.^ The quark 
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FIG. 1: The generic GIA diagram for j-3tt processes. All 
external momenta flow inward. The photon and three pions 
have momenta Q and Pi, i — 2,3, 4, respectively. 

lines in this figure represent dressed quark propagators, 
obtained as solutions of their DSE. The pion BSAs are 
solutions of the homogeneous BSE, and also the quark- 
photon vertex is dressed. 

The Lorentz structure for this diagram is 

^^(^2,^3,^4) - -ie^,p.P^PiPU{s,t,u), (6) 

where f{s,t,u) is a Lorentz scalar function of the Man- 
delstam variables s = —{Q + P2Y, t = —{Q + PzY, 
u = —{Q + Pi)"^- All three pions are on-shell, P^ = m^. 
The photon momentum is Q = —{P2 + Py, + Pi) and 
is related to the Mandelstam variables via 

s + t + u = iml-Q'^ . (7) 

The total 7 — > Stt amplitude is obtained by summing 
over all six permutations of the three pions with the ap- 
propriate charge factors. It has the same Lorentz struc- 
ture as Eq. ©, and can be written as 

Al^{P2,P5,Pi) = ^le^^p^P^pPpt F^^{s,t,u).{8) 

In the isospin limit (equal masses for up and down quarks 
and ignoring electromagnetic corrections) all six config- 
urations can be related through symmetry operations, 
which yields 

F^^[s,t,u)^ 

I (^J{s,t,u) + f{u,s,t) + f{t,u,s)) . (9) 

The process 77r — s- tttt is described by the same form 
factor: the only difference is in the kinematical variables. 
For physical three-pion production, all three Mandelstam 
variables s, t and u must be above the threshold value 
4to^. For 77r — > tttt, the kinematics are different: s > 
4 ml but t,u <0. 
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D. 



Corrections to the GIA 



Although the GIA is consistent with current conser- 
vation for three-body processes such as electromagnetic 
form factors, it is insufficient, as demonstrated in detail 
in the next section, for four-body amplitudes such as tt- 
TT scattering and F'^'^ . It was previously noted ^3 that 
if the Ball-Chiu [BC] ansatz [23 is used for the quark- 
photon vertex in combination with 



the GIA does reproduce the low-energy theorem ^ 



F^''(0,0,0) 



10.5 GeV" 



(10) 



(11) 



r^(p_l_,p) S{p) r^(p,p_). For example, for the s-channel 
gluon exchange diagrams in Fig. [21 with z = 4 and j = 3 

Go{p,P4,P3)^rAp + P4,p)S{p)T^{p,p~P3). (13) 

The corresponding amplitude for the sum of the GIA and 
the s-channel gluon exchange diagrams is (recall, momen- 
tum conservation dictates P2 — —Q — P3 — P4) 



a^{P2,P3,Pi) 



Nc [ Tr 

-^7 



S{k + Pi)G{k,Pi,P:,) 



X 5(A:~P3)r;:(fc-P3,A:-0-P3) 
xS{k-Q- Pz)T,{k -Q-P3,k + P4)] . (14) 

This treatment of the ladder diagrams is identical to the 
calculation reported for tt-tt scattering p^ . 



independent of model details for the scalar part of the 
quark self-energy, B{q^). However, this result is lost 
when the complete pion BSA (see Eq. ijJJ) is used. 

A consistent treatment requires multiple s, t and u- 
channel gluon exchange, as depicted in Fig. [3 for one of 
the six configurations. This truncation scheme is equiva- 
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FIG. 2: Gluon exchange diagrams needed to correctly de- 
scribe 7-37r processes in the rainbow-ladder truncation of the 
set of DSEs. 



III. NUMERICAL RESULTS FOR SYMMETRIC 
PION MOMENTA 

In contrast to DSE calculations for three-body pro- 
cesses, a significant computational effort is now necessary 
since G is a function of three independent variables: p^, 
p ■ Pi and p ■ Pj, for fixed Pi and Pj . Discretizing these 
variables on a three-dimensional grid, iteration is utilized 
starting with Gq. In the absence of singularities, conver- 
gence requires from 20 to 30 iterations. Close to a pole 
from an intermediate meson bound state, the number of 
required iterations increases dramatically. Typically, per 
fixed set of external variables (Q^, s, t and u) away from 
intermediate meson poles, of order 20 CPU hours are 
needed to compute the form factor to an estimated pre- 
cision of 2%. Our codes run on a parallel supercomputer 
(IBM SP) using 16 to 256 processors. Taking symmetric 
pion momenta reduces the necessary computer time by 
at least a factor of three. 



lent to summing all planar diagrams ignoring gluon self- 
interactions (triple and quartic gluon vertices). The sum- 
mation of the infinite set of ladder diagrams entails solv- 
ing an inhomogeneous BSE of the type 



Jq 



2 \ 7-^frcc 



(k) 



X V7m^('Z+) G{q, P,,Pj) S{q-)^J, , (12) 

where now q^ = q + Pi, q_ = q — Pj and 
Go {p, Pi , Pj ) is an inhomogeneous term of the form 



A. Model parameters 

The model interaction and parameters developed in 
Ref. |29J for the masses and decay constants of the light 
pseudoscalar and vector mesons are used for our numer- 
ical analysis. The effective qq interaction is 



g{k' 



47r2 D e 



■In 



1 



P/A2 



qcd; 



(15) 



^ Note that the pion decay constant depends on the current quark 
mass: the chiral hmit value is a few percent smaller than the 
physical value, = 92.4 MeV. In the model used here = 
90 MeV. 



with 7„ = 12/(33- 27V/), Nf = 4, Aqcd = 0.234 GeV, 
T = e^ - 1, T{s) = (1 - exp ■^)/s and nit = 0.5 GeV. 
The renormalization scale, /it = 19 GeV, is within the per- 
turbative domain I2J |23. The degenerate u/d quark 



mass, m 



M=lGoV 



5.5 MeV, and remaining parameters. 



uj = 0.4 GeV, D = 0.93 GeV, were determined by fitting 
"T'TTi /tt and the chiral condensate. Table Q summarizes 
TT, and predicted, p and a;, masses which are in 



fitted, , i^iiv^ j^i^v^it,uv.^, p uiivi ^j— ii" 
good agreement with observation '29]. The model's pre- 
dictions are also in very good agreement with the most 
recent JLab data [sD] for the tt^ form factor l2^. fVfQf), 
and various other pion, p and lo observables j20ll3lll33 |. 
some of which are listed in Table HI as well. In addition 
to a light pseudoscalar and vector meson, this model also 
exhibits a light scalar meson, putatively identified as the 
(7 meson [l6ll33l. 



TABLE I: Selection of calculated and measured meson masses 
and coupling constants. 





experiment [34] 
(estimates) 


calculated 
C fitted) 




5 - 10 MeV 


5.5 MeV 


- ml 


(236 MeVf 


(241 MeVt)3 




138.5 MeV 


138+ MeV 


U 


92.4 MeV 


92.5+ MeV 


mp 


771 MeV 


742 MeV 


9p 


5.01 


5.07 




6.02 


5.14 


rriui 


783 MeV 


742 MeV 


Qui 


17.06 


15.2 




400~1200 MeV 


670 MeV 



B. Quark mass dependence of 



For symmetric pion momenta, s — t ~ 
plifies to 

F3-(s,s,s) = e/(s,s,s) 



Eq. © 



(16) 



Since this reduces the necessary computer time signifi- 
cantly, we use this case to demonstrate that our scheme 
is in agreement with the low-energy theorem. 

Note that the low-energy prediction, Eq. (|ll|l . corre- 
sponds to threshold production for three massless pious 
by an on-shell photon. However, for physical pions, the 
threshold for three-pion production \s s = t = u = 4to^. 
On the other hand, for equal pion momenta, Eq. O 
reduces to = 3(m^ — s). Thus it is impossible 
to have on-shcll photons {Q^ = 0) producing three 
physical pions. Here we consider two limits, namely, 
threshold for three-pion production by timelike photons, 
F^'' {Ami, Ami, 4™^) with = _g ^2 ^ and F3^(0, 0, 0) 
with spacelike photon momenta = 3 ml . In the chi- 
ral limit, both should approach ^^^^(0,0, 0). Our ladder 
DSE result for F^'^ is presented in Fig. |21 together with 
results obtained using the GIA in combination with dif- 
ferent approximations for the vertices. There is a small 
dependence of the anomaly prediction on the current 
quark mass due to the fact that /^r depends on this mass. 
To avoid confusion with Eq. Hll|) . we denote the (mass- 



dependent) anomaly prediction by 



-^anomaly ^) 



4^2/3 



(17) 



For the physical current quark mass /^r = 92.4 MeV and 

^^a'„"omaly(0,0,0) = 9.72 GeV-^ 
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• Anomaly prediction 

o GIA, BC vertex, T = iyfi 

X GIA, BC vertex, complete pion BSA 

- + GIA, full BSE solutions for all vertices 

: * ladder DSE result 
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FIG. 3: Numerical results for the symmetric form factor 
i^^''(4m^,4m^,4m^) (top) and F^''(0,0,0) (bottom) as a 
function of the current quark mass niq divided by the model 
up/down quark mass. The curves are fits to our numerical 
results (represented by the symbols) to guide the eye. 

Our ladder DSE approach does indeed reproduce the 
correct result in the chiral limit. The strong dependence 
of the form factor at threshold, F^'^ {Ami, Ami, Ami) 
with — ~9ml, is related to the uj pole in the quark- 
photon vertex at timelike photon momentum = —m^. 
The form factor F^'^{Q, 0, 0) with spacehke Q"^ — 2i ml is 
far less dependent on the quark mass. For on-shell pho- 
tons, = 0, F^"^ {ml, ml, ml) is essentially indepen- 
dent of the current quark (and pion) mass. 

The combination of GIA, BC ansatz for the quark- 
photon vertex and a purely pseudoscalar BSA also re- 
produces the correct chiral limit. However, the functional 
behavior approaching the chiral limit value is quite dif- 
ferent from our scheme. This difference is predominantly 
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due to the lack of intermediate meson states, in partic- 
ular the which is not incorporated in the BC ansatz. 
We address the importance of intermediate meson states 
in Sec. |lV| In Fig. 01 also note that the GIA with the 
BC ansatz and full pion BSA including the pscudovector 
components F and G (see Eq. (PJ), docs not produce the 
correct chiral limit. Similarly, the GIA, with complete 
ladder BSE solutions for all vertices but without the ad- 
ditional ladder diagrams (see Fig. [5J , does not yield the 
anomaly prediction. 



C. Dependence on the photon momentum 

The dependence of the form factor on the photon mo- 
mentum is shown in Fig.^for both timelike and spacclike 
Q^, using the model quark mass corresponding to phys- 
ical pions. Again, all three on-shell pions have the same 



energy, s 



t 



— Q /3, and the threshold for 



three-pion production is Q 



^0.172 GeV^ 



25 



20 
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* ladder DSE 
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FIG. 4: Numerical results for the symmetric form factor 
F^''{s,s,s) as function of Q^. 



Our result is in good agreement with the chiral 
anomaly prediction, as is the combination of the GIA and 
the BC ansatz for the quark-photon vertex using a purely 
pseudoscalar BSA. Now the difference between the latter 
approach and our ladder DSE result is more apparent 
as our result displays a much stronger Q^-dependence. 
In particular, our result indicates a pole-like singularity 
for timelike Q^, which can be attributed to the presence 
of a vector meson pole in the dressed quark-photon ver- 
tex. This singularity is evident in both the GIA and 
our full calculation, but not in the calculation using the 
BC ansatz. Again, the GIA does not reproduce the 
anomaly prediction without the additional ladder dia- 
grams of Fig. [3 



IV. COMPARISON TO MESON EXCHANGE 
MODELS 

The ladder DSE result clearly indicates a pole-like sin- 
gularity in the timelike region, due to the presence 
of a vector meson pole in the dressed quark-photon ver- 
tex. In addition to this pole, there are pole-like singu- 
larities from intermediate meson exchanges in the s, t 
and u channels, similar to the scalar and vector meson 
exchange contributions in tt-tt scattering Here how- 
ever, the scalar meson exchange contribution is forbidden 
by spin/parity conservation. Thus only vector meson in- 
termediate states contribute to 7-37r processes. 



A. Vector meson dominance 

The effective vector meson dominance [VMD] La- 
grangian is 



(18) 



9p 9u. 

where tt, /v, 'j-'i/ and are the tt, p, lu and photon fields, 
muj, mp are the uj, p masses and gujp-w, Qp-n-K are the w 
pn, p — + TTTT coupling constants respectively. With this 
Lagrangian the VMD form factor for 7 ^ cj ^ pTr — + Stt, 
is I35II 



VMD 



2 6771^ QujpTTQpTTTT 



1 



— t 



(19) 



Note that this VMD form factor has vector meson poles 
in both the photon momentum and the Mandelstam 
variables s, t and u. The pole in is associated with 
an intermediate w meson, whereas the poles in s, t and 
u are associated with p intermediate states. 
The measured uj, p e^e~ decay widths 



ee 



Ana rriu, 

^gT 



0.60 keV, 



6.85 keV, 



(20) 
(21) 



determine the w-photon, g^j — 17.06, and p-photon, 
gp = 5.01, couplings which compare reasonably with the 
model calculations of 15.2 and 5.07, respectively 0| (see 
also Table nj. The hadronic coupling, gp-^Tn was previ- 
ously predicted ^32] in this model to be 5.14 which agrees 
favorably with the experimentally extracted value of 6.02. 

Due to limited phase space, the lo pi: decay is forbid- 
den, precluding direct determination of gi^pT^- However, 



this coupling can be related to the uj 
using VMD for u; 



ag, 



2 {ml 



7r°7 decay, again 



ml) 



2 N3 



(22) 



6 



Using the measured width, F"^ = .734 MeV, yields 

= 11.8 GeV"^ With these parameters the VMD 
prediction for the 7r° — > 77 anomaly, via tt" up 77, 



9i 



.0253 GeV" 



(23) 



This is in excellent agreement with both the low-energy 
theorem, ^2-^(0) = a/inf^) = .0251GeV"\ and the ex- 
perimental value, .025 ± 0.001 GeV"^ 

Returning to the F^"^ {s,t,u) form factor, the VMD 
prediction in the chiral and zero momenta limits is 



B. Comparison of DSE and VMD for 7 



Because of the established phenomenological success 
of VMD, it is worth while to compare the two VMD 
form factors to our DSE result. In particular, it is of 
interest to ascertain that the DSE approach can repro- 
duce the resonance features of the VMD. In doing so, 
consistency mandates that the VMD phenomenological 
parameters be replaced by the corresponding values cal- 
culated in the DSE model. All necessary DSE coupling 
constants, except for g^jpir which is discussed below, are 
listed in Tabic Q] In Fig. |Slthc pure, Eq. lfTO|l . and mod- 



VMD 



(0,0,0) 



6 G Qujp-KQ p-K1T 



12.6 GcV" 



(24) 



This is above the low-energy theorem result, Eq. (|ll|l . 
but in potentially better agreement with the observed 
value, 12.9 ± 0.9 ± 0.5 GeV~'^, measured at higher energy 
(note that Ref. [s^ quotes a corrected experimental value 
of 11.9 GeV~^). By reducing gp^^ it is possible for the 
VMD form factor to agree with the low-energy theorem, 
however this conflicts with the KSRF relation j37t i38| | 



9p7rir = mp/{V2f.„). 



(25) 



Further, the VMD result is not consistent with anoma- 
lous Ward identities of Ref. ^ . These two issues can be 
resolved by adding a contact term to the Lagrangian 



G^e^"^"^,.V,7r • Va^ x V^t? , (26) 



with the coupling G^^ governing the direct uj ^ Sir decay. 
This gives the modified VMD result pHiof 



—6 e m 



9u. 



+ mf 



3tt 

VMD 



{s,t,u) . 



(27) 



References j^E 1^ determine Gu, by satisfying both 
Eq. H25|l and the low-energy theorem for 7-37r. Further, 
these analyses also invoke SU(2) flavor symmetry and 
universality using the relations g^j = Sgp = SgpT^r, which 
are reasonable, to obtain 



Gf, 



9 PTTTT 

16^ 



(28) 



With these relations the modified form factor can be 
more succinctly expressed in terms of only /tt, using 
Eq. l(T7|) , instead of meson coupling constants 



^2 ^2 ^ 



t 



(29) 
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ladder DSE 

pure VMD, fit near pole 
modified VMD 
QHD fit 



photon Q [GeVl 

FIG. 5: Comparison of the ladder DSE calculation (aster- 
isk) to the DSE parameterized pure (dotted) and modified 
(dashed) VMD form factors, F^'^ (s, s, s), as function of Q^. 
The solid curve represents a phenomenological QHD fit to the 
DSE form factor. 

ified, Eq. (|29|l . VMD form factors are compared to the 
DSE result from « to the uj resonance region. 

The ladder DSE and VMD results are in good qualita- 
tive agreement. In particular, the lu resonance exhibited 
in the VMD is well reproduced by the ladder DSE ap- 
proach. Note that the DSE approach naturally agrees 
with the low-energy theorem while only by construction 
can the modified VMD form factor reproduce this result; 
the pure VMD form factor does not reproduce the correct 
low-energy limit. 

Because of off-shell ambiguities, one can not consis- 
tently calculate the unphysical Lopir amplitude. However, 
by fitting the pure VMD form factor, Eq. H19() . to the 
DSE result near the uj pole and using the calculated val- 
ues of Table for the remaining coupling constants in 
that equation, we extract the value g^jpir — 10.3GeV"\ 
This is slightly below the phenomenological VMD result 
of 11.8 GeV^^ but significantly lower than other theoret- 
ical values: SU(3) flavor symmetry, IBGeV"^ [ill; QCD 
sum rules, 15 to 17GeV~^ 0; hght-cone sum rules, 
IGGeV"^ m. 

In addition. Fig. |S1 also displays a form factor based on 
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a purely phenomenological quark-hadron model [QHD] 
incorporating both contact interactions and meson- 
exchange contributions 



6 m^, 



'37T 



5c. (Q2 + ml) 



-ynp 



(30) 



The three coefficients C are determined by fitting the 
ladder DSE calculation at different s, t, u and values. 
An accurate fit over a large kinematical range is obtained 
with C^3^ 0.248 GeV-^ C^a^r = 5.94 GeV"^ and 
C-yi^p = 0.028 GeV~^. This provides a useful analytic 
representation of the DSE form factor. 



C. Comparison of DSE and VMD for 77r -kit 

Finally, we consider the process 77r — > tttt. It is de- 
scribed by the same form factor F^^ {s,t,u) as the pro- 
cess 7 Stt, but the kinematical variables are different. 
In general, s ^ t ^ u, but for simplicity, we only con- 
sider the more symmetric case t ~ u. The kincmatically 
allowed region is s > 4to^ and t,u <Q. 
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FIG. 6: Comparison of the ladder DSE (asterisk), pure VMD 
(dotted), modified VMD (dashed) and fitted QHD (solid) 
form factors, F'^^{s,t,t), as function of s at = 0. 

In Fig. El we compare our DSE results with the VMD 
and phenomenological meson exchange form factors. For 
fixed = 0, the uj meson does not appear as a pole in 
the form factor. The p pole in the s channel on the other 
hand is clearly identifiable. In addition, there are also 
poles in the unphysical region, corresponding to interme- 
diate p mesons in the t- and w-channels. 

With the same parameters as in the previous subsec- 
tion, the pure VMD form factor agrees very well with the 



ladder DSE calculation near the p pole, but overestimates 
the form factor by about 15% near threshold. The modi- 
fied VMD form factor on the other hand agrees very well 
with the DSE result near threshold, but starts to deviate 
for increasing values of s. The purely phenomenologi- 
cal QHD form factor can describe the DSE form factor 
accurately both near threshold and in the resonance re- 
gion, not only at = but also at other values of Q^. 
However, this is not surprising, since it was fitted to our 
DSE result to provide an analytic representation of our 
results. 



COMPARISON WITH EXPERIMENTAL 
DATA 



Finally, for the benefit of experimentalist the same 
three form factors are compared to the limited data. We 
use the phenomenological parameters in the VMD form 
factors permitting their most favorable prediction. In 
particular, the p mass in Eqs. (I19|l and H29|) is replaced 
by nip — IFp where Fp = 149 MeV is the experimental 
width of the p meson. This changes the p-meson pole 
singularity in Fig. El to a resonance peak with a finite 
maximum. We ignore the much smaller lu width since 
we consider the process 77r —f tttt for on-shell photons 
(g2 = 0) only. 

The ladder truncation of the inhomogeneous BSE 
Eq. (|12|l . on the other hand, has a real pole singularity at 
the p-mass. One would need to include two-pion inter- 
mediate stated in the interaction kernel to consistently 
produce a p-meson width. 
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FIG. 7: Comparison of DSE (asterisk), pure VMD (dotted), 
and modified VMD (dashed) form factors, F^^ {s,t,t), with 
data at — 0. The experimental coupling constants, masses 
and width are used in both VMD form factors (the shaded 
region is unphysical) . The data point is from Ref. (re- 
analysed in Ref. pal): new experimental results from JLab 
are anticipated in the region 0.27 GeV^ < s < 0.72 GeV^ 
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In Fig. [7| we display the absolute value of both the 
pure and the modified VMD form factors as functions 
of s, the total energy of the two outgoing pions, for t = 
u = (3m^ — s)/2 and on-shell photons. For comparison, 
we also show the ladder DSE calculation below the p- 
pole singularity. The modified VMD form factor appears 
to be the most realistic and agrees quite well with the 
microscopic DSE calculation near threshold, s = 4m^. 
At larger values of s it reflects the p-meson width and 
therefore exhibits a resonance peak rather than the pole 
singularity of the ladder DSE calculation. The large error 
bar on the experimental data point precludes definitive 
conclusions but new, more precise measurements are in 
progress which should provide useful insight and model 
constraints. 



VI. CONCLUSION 

In summary, the anomalous form factor for "f-Sn pro- 
cesses has been calculated in a self-consistent Dyson- 
Schwinger, Bethe-Salpeter formulation in the rainbow- 
ladder truncation. The significant finding is that to re- 
produce the fundamental low-energy theorem, it is nec- 
essary to go beyond the generalized impulse approxima- 
tion. Specifically, one has to include complete sets of 
ladder diagrams in the s, t and u channels. 

The ladder DSE result is in good agreement with vec- 
tor meson dominance. The ladder dressing of the quark- 
photon vertex is essential for generating the w resonance 



in the process 7 — > w — > Stt. Inclusion of the ladder 
diagrams beyond the GIA (see Fig. ^ is necessary to 
produce an intermediate p-meson state in the two-pion 
channel. A meson-exchange form factor incorporating 
both contact interactions and meson-exchange contribu- 
tions, fitted to our DSE calculation, provides a useful 
analytic representation of our results over a large kine- 
matical range. 

Our predictions are also in agreement with the lim- 
ited data near threshold and await confrontation with 
additional and more precise measurements currently in 
progress. For larger values of s we expect the modified 
VMD form factor to be in better agreement with exper- 
iment since it includes the /9-meson width. These results 
should also be of interest to experimentalists investigat- 
ing double pion photo-production. 

Future work will address calculating the width for 
u —^ Sir as well as the strangeness processes 7 KKt: 
and "fK ttK. Longer term, this framework will be ex- 
tended to describe meson electro-production from a nu- 
cleon target. 
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